Spectroscopy of a cold strontium Rydberg gas 
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We present a study of a cold strontium Rydberg gas. The narrowband laser excitation of Rydberg 
states in the range n = 20 — 80 from a 6 mK cloud of strontium atoms is detected using the 
spontaneous ionization of the Rydberg atoms. Using a high-resolution step-scanning technique, we 
perform detailed measurements of the Stark maps of selected Rydberg states. We find excellent 
agreement between the measured Stark maps and a numerical calculation based on an independent- 
electron model. Finally we show that excitation of the second valence electron can be used to probe 
the dynamics of the Rydberg gas with nanosecond temporal resolution via autoionization. 



INTRODUCTION 

Rydberg atoms offer an ideal opportunity to study 
strongly interacting quantum systems, due to the long- 
range dipole-dipole interactions between them. Two in- 
dividually trapped neutral atoms excited to Rydberg 
states have been shown to exhibit a high degree of en- 
tanglement, through the dipole blockade mechanism 0- 
4]. These entangled states are important for the study 
of quantum information, as they can be used to per- 
form quantum gate operations [5J. In mesoscopic atom 
clouds, the dipole blockade effect leads to the formation 
of highly-correlated many-body states @ that can exhibit 
long-range crystalline order 0, Q . The dipole blockade 
also modifies the atom-light interactions Q , giving rise to 
cooperative optical effects which could be exploited 
to create photonic phase gates [ll|. The long-range in- 
teractions cause Rydberg gases to spontaneously evolve 
into cold neutral plasmas |12j | through spontaneous ion- 
ization. 

Most experiments in this field have used alkali metal 
atoms, where only a single valence electron is available for 
manipulation. Divalent atoms, such as strontium, open 
up new opportunities for the study of Rydberg gases. For 
Rydberg states of low angular momentum L, excitation of 
the inner valence electron leads to autoionization, which 
can be used as a state-selective probe of the interactions 
in a cold Rydberg gas Off resonance, the strong 

optical transitions of the ionic core can be exploited to 
create optical dipole traps for Rydberg atoms [l4j. For 
states of high L the inner valence electron can be excited 
without ionizing the atom [l5j . raising the possibility of 
imaging the Rydberg atoms, in a similar manner to the 
direct imaging of Sr + ions in a cold plasma [l6j . 

In this paper we discuss the creation and manipulation 
of a cold strontium Rydberg gas. We describe the cold 
Rydberg gas apparatus (section [J) , and present a high- 
resolution spectroscopic study of selected Rydberg states 
at principal quantum number n > 50 (section [TTJ) . In or- 
der to understand the behaviour of cold Rydberg gases, 
a clear understanding of the nature of the dipole-dipole 
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FIG. 1. (a) A scale diagram of the vacuum chamber. Zeeman- 
slowed atoms enter from the left. Electrodes create an electric 
field E that directs ions to the MCP. The fluorescence from 
the MOT is monitored by a photodiode (PD). The Rydberg 
excitation beams Ai and A2 (see figure [2j a)) are also shown, 
(b) A schematic of the apparatus inside the vacuum cham- 
ber. The MOT coils are wound on copper formers, which also 
support eight electrodes in a split-ring geometry. 



interactions is essential. This requires a knowledge of the 
dipole matrix elements. We present a detailed numerical 
calculation of the Stark maps based on dipole matrix el- 
ements calculated using an independent electron model, 
and find very good quantitative agreement with experi- 
ment (section [TTT|) . Finally we show that laser excitation 
of the second valence electron can be used as a high time- 
resolution probe of the Rydberg gas (section IIV|) 



I. EXPERIMENTAL SET-UP 

Cold 88 Sr atoms are prepared in a standard six-beam 
magneto-optical trap (MOT) operating on the 5s 2 1 So — > 
5s5p 1 Pi transition at 461 nm. This light is produced by a 
frequency doubled diode laser system (Toptica TA-SHG). 
To stabilize the laser at 461 nm, we perform polarization 
spectroscopy [l7j in a strontium dispenser cell [l8j . A 
schematic of the trapping region is shown in figure [lja). 
A MOT for strontium requires a magnetic field gradi- 
ent of ~ 30 Gem -1 . The field coils are mounted inside 
the vacuum system, so that this field can be produced 
without the need for water cooling. The coils are made 
from Kapton insulated copper wire wound onto copper 
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formers. 

A micro-channel plate (MCP) detector is mounted in- 
side the vacuum system to detect ions produced in the 
Rydberg gas. The MCP is positioned away from the line 
of-sight of the atomic beam, as shown on figure HJa), 
and is further protected by a baffle. A wire mesh grid 
electrode placed in front of the MCP is used to control 
the stray electric field from the detector. The grid is 
typically held at —10 V. An additional electric field is 
used to direct ions to the MCP. The field is created by 
eight independently controllable electrodes, arranged in 
a split- ring geometry The electrodes are mounted 

onto the MOT coil formers, as shown in figure [IJb), and 
insulated from the copper with ceramic spacers. 

The trap is loaded from a Zeeman slowed atomic beam. 
Strontium is heated in an oven, which consists of heater 
wire wrapped around the outside of a standard stainless 
steel vacuum tube [2(J ■ The hot strontium is collimated 
using a bundle of 8 mm long, 170 /jm diameter steel cap- 
illaries, held in a nozzle which is separately heated to 
prevent condensation of the strontium. The hot atomic 
beam passes through a 30 cm long solenoid-type Zeeman 
slower. The Zeeman slower is encased in a mild steel yoke 
to boost the magnetic field at the ends of the slower, and 
further shield the trapped atoms from stray magnetic 
fields. We typically trap 4 x 10 6 atoms, at a density of 
2xl0 10 cm~ 3 and a temperature of 6 mK. 

The Rydberg sample is prepared using a two-step ex- 
citation Ai + A2, as illustrated in figure^ a). The light 
for the first step, Ai = 461 nm, is derived from the cool- 
ing laser, and is resonant with the 5s 21 So — > 5s5p 1 P\ 
transition. The light for the second step (the Rydberg 
excitation), A2, is produced by one of two frequency dou- 
bled diode laser systems (Toptica DLS SHG), at 420 nm 
or 413 nm, depending on the range of Rydberg states 
we wish to access. Both beams are pulsed on simulta- 
neously for 4 /is. The beams are overlapped, counter- 
propagating, and linearly polarized in the vertical direc- 
tion. They propagate normal to the direction of the elec- 
tric field at the position of the cold atoms. 

The wavelength of the fundamental of the Rydberg 
excitation laser is measured using a High Finesse WS7 
wavemeter, and is converted into the output laser fre- 
quency U>2- 

The next section outlines a method for performing Ry- 
dberg state spectroscopy using this experimental set-up. 



II. HIGH RESOLUTION RYDBERG 
SPECTROSCOPY 

We perform spectroscopy by varying the frequency 
L02 of the Rydberg excitation laser, and measuring the 
amount of spontaneous ionization. Spontaneous ioniza- 
tion has been observed in nearly all experiments on cold 
Rydberg gases. The dominant mechanism is thought to 
be interaction-enhanced Penning ionization [2l], [22[ • In 
our experiment ionization could also arise from collisions 
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FIG. 2. (a) Energy level diagram for the two-step excitation, 
Ai + A2, to a Rydberg state 5snL (b) Timing diagram for 
the excitation sequence. The Rydberg state is populated af- 
ter the ground-state atoms are released from the MOF. The 
MOT light is pulsed before and after the Rydberg excitation, 
(c) Experimental sequence for the step-scan technique. The 
Rydberg excitation sequence is repeated 100 times. The aver- 
aged ion signal and fluorescence, and the frequency u>2 of A2 , 
are recorded. The frequency U12 is stepped, and the process 
repeated to produce a spectrum. 



with hot atoms not captured by the Zeeman slower [12j . 
We observe ionization for all of the states we have stud- 
ied, which include states of principal quantum number 
n = 18 and above. The ionization is present regardless 
of whether the interatomic interactions are attractive (as 
is the case for the 1 So states in Sr [lU), or repulsive (as is 
the case for the X D2 states above n = 25). The ionization 
is rapid |23J ; we observe ionization within 1 /xs of starting 
the Rydberg excitation. At the Rydberg densities used in 
this paper (< 5 x 10 s cm -3 ), the spontaneous ionization 
signal is proportional to the Rydberg excitation laser in- 
tensity. At higher Rydberg densities we observe the for- 
mation of a cold plasma 1131 . which leads to avalanche 
ionization 
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The spontaneously created ions are directed to the 
MCP with a4/is long electric field pulse. The amplitude 
of this pulse (4 Vcm -1 ) is not sufficient to field ionize 
any of the Rydberg states described in this paper. The 
signal from the MCP is amplified, recorded on a digital 
oscilloscope and integrated to yield a total ion signal in 
V/iS. 

The sequence that we use to perform spectroscopy is 
illustrated in figure EJb) . A pulse of MOT light before 
and after the Rydberg excitation is used to measure the 
ground state atom number and any loss of atoms. This 
excitation sequence is repeated 100 times, and the ion- 
ization and fluorescence signals are averaged on the os- 
cilloscope. After the 100 repeats the average signals are 
recorded, and the frequency W2 is sampled 100 times. A 
voltage is applied to the laser piezo to step the frequency 
ll>2, and whole process is repeated, as illustrated in fig- 
ure [DJc). By repeatedly stepping u>2 we can build up a 
Rydberg state spectrum. We refer to this method as the 
"step-scan" technique. The entire stepping and recording 
process is automated under computer control. 

Importantly, measuring u>2 at each step enables us to 
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FIG. 3. Example spectra obtained using the step-scan tech- 
nique. The detunings Ai, A2 are defined relative to our mea- 
sured line centres, (a) Saturated absorption spectrum of 
the 5s 2 1 So — > SsSp 1 ?! transition. The fit is a sum of six 
Lorentzians. (b) Spectrum of the 5sl9d 1 D2 state. The solid 
line is a Gaussian fit with a FWHM of 48 MHz, which matches 
the linewidth of the 5s5p 1 Pi state (32 MHz) with expected 
power-broadening. 
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FIG. 4. Spectrum in the region n = 80, in the presence of 
a 0.35 Vcm -1 electric field. The detuning A2 is relative to 
the energy of the 80 1 T>2 state in zero field. The 80 1 T>2 state 
signal peaks at 1.4 V/xs. The electric field allows otherwise 
forbidden transitions, and has split the 80^1 and 80 3 D2 
states into their \mj\ components. The lines are assigned 
using the single electron simulation. 



negate the effect of slow drifts in the laser frequency. 
To characterize the accuracy of this step-scan technique, 
we performed saturated absorption spectroscopy of the 
5s 2 x So — ► SsSp 1 ?! transition, as shown in figure 
The frequency u>% of the 461 nm laser is stepped across 
the resonance in steps of 1.4 MHz, which are clearly re- 
solved by the wavemeter. We fit the spectrum with a 
sum of six Lorentzians, one for each isotope and hyper- 
fine component [24]. The free parameters are the overall 
amplitude, a sloping background and a frequency axis 
scaling parameter a, such that lo\ — > au}\. From the fit, 
we obtain a = 1.03 ± 0.02, indicating that the relative 
frequency accuracy is very good. By comparing our mea- 
sured wavelengths to literature values, we estimate that 
the absolute accuracy of the wavemeter is ~ 200 MHz. 
This is consistent with the wavemeter specification when 
using multi-mode optical fibres. 

Examples of the Rydberg spectra that we obtain are 
shown in figures E^b) and 2] Figure Elb) shows a high- 
resolution scan across a single Rydberg state. The 
linewidth of the Rydberg excitation is consistent with the 
power-broadened linewidth of the intermediate SsSp 1 ?! 
state. In figure Q] the Rydberg laser is stepped over 
8 GHz, in steps of 7 MHz, in the vicinity of the 5s80d 1 T>2 
state. This spectrum was acquired in the presence of 
a constant electric field of 0.35 Vcm -1 , which allows 
the excitation of the dipole-forbidden 1 Pi state. Triplet 
states are also visible in the spectrum, indicating that 
singlet-triplet mixing [25| is still present even at n = 80. 
The 81 3 Si peak is ~ 200 times smaller than the principal 
80 1 D2 peak, which highlights that spontaneous ioniza- 
tion is detectable even for these extremely weakly popu- 
lated states. 



III. SINGLE ELECTRON MODEL AND STARK 
MAPS 

Rydberg atoms exhibit greatly enhanced interatomic 
interactions as compared to ground state atoms. To cal- 
culate the size and form of these interactions the dipole 
matrix elements between different atomic states must be 
known. This requires a knowledge of the atomic state 
wavefunctions. For the alkali elements, the waveftmc- 
tions are commonly calculated by considering a single 
electron moving in a model potential, which incorporates 
the effect of the valence electron penetrating the closed 
electronic shells [26j . 

The situation for alkali earth elements is more com- 
plex, since there are two valence electrons. This leads 
to perturbations in the energy level structure due to the 
presence of doubly excited perturber states [13] ■ To fully 
describe the energy level structure of the alkali earth el- 
ements, multi-channel quantum defect theory (MQDT) 
must be used [28| . However, [29| recently obtained good 
agreement between experimental and theoretical Sr Stark 
maps in the region of principal quantum number n = 12, 
using dipole matrix elements obtained by treating the 
two electrons as completely independent and moving in 
a simple model potential. In the following, we extend this 
model to much higher values of n and total orbital angu- 
lar momentum L, and find that it remains in good agree- 
ment with our experimentally measured Stark maps. 

Our simulation uses known quantum defects [2o1 . 130L 
l3lj . and extrapolates them to regions where they have 
not been previously measured. Using these values we fit 
a model potential of the Klapisch form [32j using a simu- 
lated annealing algorithm. Radial wavefunctions are cal- 
culated, using the Numerov method [33], by considering 
a single electron moving in this model potential. Dipole 
matrix elements are calculated from the radial waveftmc- 



4 




E (Vcrrf 1 ) 

FIG. 5. Stark map in the region of the 5s56d 1 T>2 state. The 
detuning A2 is relative to the energy of the state in zero elec- 
tric field. The blue (solid), red (dashed), green (dot-dashed) 
lines are the simulated \mj\ = 0, 1, 2 components respectively. 
The grey region is the high angular momentum manifold. The 
inset shows the fit of the simulated Stark shift averaged over 
the three \m.j\ components (red line) to the experimental cen- 
tre frequency (black dots), extracted using a Gaussian fit to 
the lineshape at each electric field. The linewidth (FWHM) 
is indicated by the error bars. 

tions [34| . To calculate a Stark map in the region n = no 
states of n — tiq ± 5 are included, and the angular mo- 
mentum L states are truncated at L — 40. 

The Stark map is measured using the step-scan tech- 
nique, with a static electric field applied during the Ry- 
dberg excitation. Measured and simulated Stark maps 
are shown in the region of the 5s56d 1 D2 state in fig- 
ure 03 and the 5s80d X D2 in figure |5] In order to compare 
experiment and theory, the calibration factor between 
the applied voltage and the electric field must be deter- 
mined. The calibration factor is obtained by fitting the 
simulated Stark map to the experimental data, as shown 
in the inset to figure [5j There are three fitting param- 
eters: the calibration factor, the voltage offset due to 
the stray electric field, and a frequency offset due to the 
wavemeter uncertainty. We constrain the voltage offset 
using the variation in the amplitude of the 80 : Pi peak 
with electric field. This state cannot be populated in 
zero field, since we excite via the SsSp 1 ?! state. The fit 
shown in the inset to figure [5] yields a calibration factor 
of 0.205 ± 0.005 Vcm -1 per volt applied. A similar fit 
to the n = 80 Stark map yields a calibration factor of 
0.215 ± 0.005 Vcm -1 per volt applied. The correspond- 
ing value for the stray electric field obtained from the 
voltage offset is 2.8 ± 0.4 mVcm -1 . This is consistent 
with the estimated electric field due to the MCP grid. 

These calibration curves highlight the excellent quan- 
titative agreement between the single electron model 
and the experiment for the 1 D2 states. In figure [SJa), 
the splitting between the \m,j\ — 0,1,2 (blue solid, red 
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FIG. 6. A Stark map in the region of the 5s80d 1 T>2 state. 
Solid black lines are the experimental spectra. The detuning 
A2 is relative to the energy of the 80 X D2 state in zero electric 
field. The blue (solid), red (dashed), green (dot-dashed) lines 
are the simulated \mj\ — 0, 1, 2 components respectively, and 
the grey shaded areas represent the simulated high angular 
momentum manifolds. The 80 3 D3 state is doubly forbidden 
and is not observed. 



dashed, green dot-dashed lines) agree well with the data 
(black lines) , and the interaction with the higher angular 
momentum manifold (grey region) is correctly predicted 
and visible in the data. The agreement is also very good 
across the n — 80 Stark map shown in figure H3 The only 
previously measured state energy in this region is for the 
5s80d 1 D 2 state [3(J, and we have adjusted the quantum 
defects for the other series to match our measured split- 
ting from the 80 X D 2 state in zero field. This requires a 
change of < 1% from the extrapolated quantum defects. 
The interaction of the low L states with the high angu- 
lar momentum manifolds (shaded grey areas on figure 
is clearly visible. The agreement is better for the sin- 
glet states than the triplet states. More quantum defects 
must be found by extrapolation for the triplet series, as 
the spectroscopic data is less complete. Alternatively, the 
disagreement may reflect the fact that the single electron 
model ignores the mixing between the singlet and triplet 
states. 

We have shown that a single electron model can be 
used to calculate dipole matrix elements for strontium, 
over a large range of n. The single electron model is 
therefore appropriate for calculating the dipole-dipole in- 
teractions, which depend on these matrix elements flij . 
The single electron model does not work for predicting 
state lifetimes and line strengths, which are strongly per- 
turbed by the presence of doubly excited states. To fully 
describe the energy level structure of the alkali earth el- 
ements multi-channel quantum defect theory (MQDT) 
must be used [28j . which takes into account the divalent 
nature of the atom. 
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FIG. 7. The autoionization signal for the 5sl9d D2 state as 
a function of time. At = is defined as the end of the 4 fis 
long excitation to the Rydberg state Ai + A2, as illustrated in 
the inset. The solid line is the solution to the optical Bloch 
equations. 
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FIG. 8. Autoionization spectra of the (a) 5sl9d D2 and (b) 
5s20s 1 So states. The spectra are fit with a six-channel MQDT 
model. A3 is the detuning of W3 from the 5si/2 5p3/2 ion 
transition. These spectra are taken at a delay At = 0.5 lis. 



IV. PULSED AUTOIONIZATION 

Rather than relying on spontaneous ionization, we can 
detect the Rydberg atoms using autoionization by ex- 
citing the inner valence electron. If an atom is initially 
prepared in a Rydberg state, the inner valence electron 
can be excited independently, a process known as iso- 
lated core excitation . The doubly excited states cou- 
ple not only to each other (as with singly excited states), 
but also to the degenerate continuum. For states of low 
angular momentum L this leads to rapid ionization, and 
these states are referred to as "autoionizing" . There has 
been much research into the nature of autoionizing res- 
onances, with MQDT proving to be extremely powerful 
for describing doubly excited states [28| . The autoioniza- 
tion process has a large cross-section, and hence produces 
significantly more ions than spontaneous ionization. We 
typically find that by exciting an autoionizing resonance 
the signal is increased by a factor of 100. 

We prepare atoms in the 5snZ Rydberg state, using 
the same sequence as before (figure [5]) . At a variable 
time At after the end of the excitation to the Rydberg 
state, a 10 ns light pulse at A3 = 408 nm excites the inner 
valence electron, creating the autoionizing state 5p 3 / 2 n'Z. 
The light A3 is provided by a tunable pulsed dye laser 
system, with a bandwidth of a few GHz, and a repetition 
rate of 10 Hz. The autoionizing pulse is immediately 
followed by an electric field pulse that directs the ions to 
the MCP. 

At each value of the wavelength A3 and delay At we 
measure a Rydberg spectrum using the step-scan tech- 
nique described previously The spectrum is normalized 
by the average pulse intensity measured on a photodiodc, 
and the ground state atom number. From this spectrum 
we extract the peak autoionization-enhanced ion signal, 
S, in units of V/zs/iW -1 cm 2 per 10 6 atoms. 

By fixing the wavelength A3 and scanning the delay At 
across the 4 /zs Rydberg excitation pulse, the evolution of 



the Rydberg population can be mapped with high time 
resolution as shown in figure [JJ Both the Rydberg excita- 
tion and subsequent decay are clearly visible. By fitting 
an exponential to the decay we find a state lifetime of 
7"i9D = 880 ± 30 ns. The complete time evolution can be 
modelled by solving the three level optical Bloch equa- 
tions for the scheme 5s 21 So -> 5s5p 1 Pi — > 5sl9d 1 D 2 . 
The Rabi frequencies and decay constants for each step 
of the excitation are constrained using our measurements 
of the laser intensity and the lifetime T19D, and the os- 
cillator strength for the 5s5p 1 Pi — > 5sl9d 1 D 2 transition 
from [36|. The only remaining free parameter is the am- 
plitude scale. The fit is in very good agreement with 
the data, as can be seen in figure [7J This indicates that 
the autoionization signal is proportional to the Rydberg 
state population. An identical analysis for the 5s20s 1 So 
state yields a state lifetime of 3.1±0.1 /is. 

By fixing the delay At, and varying the wavelength A3, 
the spectrum of the 5snl — > 5p 3 / 2 n'Z autoionizing transi- 
tion can be measured. MQDT allows us to reproduce the 
spectra of the autoionizing states, which in turn allows 
us to identify exactly which state a sample is in from its 
autoionization spectrum. The spectra are fit with a six- 
channel MQDT model, using the same fitting parameters 
as [37[ for the 19 1 D 2 state (with the ellip tical polariza- 
tion of A3 included in the model), and [38[ for the 20 1 So 
state. The difference in the shape of the spectra is due 
to the difference e = <5& — S a between the bound (Ryd- 
berg) state quantum defect 5b, and the autoionizing state 
quantum defect 5 a . For e ~ (modulo-one) each Ryd- 
berg state only overlaps with a single autoionizing state, 
and the autoionizing spectrum is a single peak, as seen 
for the 20 % state with e = 0.15. For e ~ 0.5 (modulo- 
one) each Rydberg state overlaps with two autoionizing 
states, and the autoionizing spectrum is double-peaked, 
as seen for the 19 1 D 2 state with e = 0.46 

The spectra in figure [5] illustrate that autoionization 
can be used for state-selective detection. It is particu- 
larly sensitive to changes in L. In previous work we have 
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identified, and quantified, transfer of population from one 
Rydberg state to another through analysis of the autoion- 
ization spectra fl3Tj . In the region n « 20 the autoion- 
ization spectra are approximately 1 THz wide as shown 
in figure [8] Even at n — 56, the width of the autoion- 
izing spectrum for low-L states is > 10 GHz 13]. The 
high temporal resolution of the autoionization technique 
displayed in figure[7]can therefore be fully exploited with- 
out compromising its state-selectivity. The combination 
of high temporal resolution, state-selectivity and large 
signal-to-noise ratio make autoionization a powerful tech- 
nique for probing the dynamics of cold Rydberg gases. 

CONCLUSION 

The use of alkali earth elements such as strontium of- 
fer new ways to study and manipulate Rydberg atoms. 
We have performed the first study of Rydberg states in a 



cold gas of strontium. We have used, and experimentally 
verified, a simple single-electron model to find dipole ma- 
trix elements, which can be used to calculate interatomic 
interaction strengths. We have also used autoionization 
through two-electron excitation as a high-yield probe of 
a cold Rydberg gas, which can measure state population 
dynamics with high temporal resolution. By focusing 
the autoionizing laser this technique could be used to 
spatially probe the Rydberg gas, on length scales com- 
parable to the dipole blockade radius. 

We thank J.D. Pritchard for assistance with the dipole 
matrix element calculations, and C. S. Adams and D. 
Carty for the loan of equipment. This work was sup- 
ported by EPSRC Grants No. EP/D070287/1 and No. 
LP/82000, and by Durham University. 
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